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I.  STATEMENT  OF  THE  PROBLEM 

The  propagation  of  electromagnetic  waves  in  snow  at  microwave  and 
millimeter  wave  frequencies  is  affected  by  volume  scattering.  Snow  is  a 
dense  nontenuous  medium,  in  a  dense  medium,  the  particles  occupy  an 
appreciable  fractional  volume.  In  a  nontenuous  medium,  the  permittivities 
of  the  particles  are  significantly  different  from  that  of  the  background 
medium.  Our  studies  are  directed  towards  basic  theoretical  investigations 
on  effects  of  multiple  scattering  of  waves  in  dense  nontenuous  media.  We 
apply  the  theory  to  study  the  characteristics  of  coherent  and  Incoherent 
waves  in  snow  at  microwave  and  millimeter  wave  frequencies. 


II.  SUMMARY  OF  IMPORTANT  RESULTS 

A.  Dense  Media  Radiative  Transfer  Equation 

Snow  can  be  characterized  as  dense  media.  We  have  made  progress 
in  the  study  of  active  and  passive  microwave  remote  sensing  of  dense  media. 
In  a  dense  nontenuous  medium,  the  classical  radiative  transfer  theory, 
which  is  based  on  the  assumption  of  independent  scattering,  is  not  valid. 

A  set  of  transfer  equations  has  been  derived  for  dense  nontenuous  media. 
The  derivation  is  based  on  field  theory  under  the  quasi-crystalline 
approximation  with  coherent  potential  on  the  first  moment  of  the  field  and 
the  modified  ladder  approximation  on  the  second  moment  of  the  field.  These 
equations  are  called  radiative  wave  equations.  They  Include  a  summation  of 
all  the  ladder  terms  and  assume  a  form  that  is  identical  to  the  classical 
transfer  equations.  However,  the  relations  of  the  extinction  coefficient, 
the  scattering  coefficient,  the  albedo,  and  the  phase  functions  to  the 
physical  parameters  of  the  medium  are  modified  to  include  the  effects  of 
dense  media.  The  new  equations  still  preserve  the  advantages  that  (l) 
multiple  scattering  of  the  incoherent  intensities  are  included,  (ii)  energy 
conservation  and  reciprocity  are  obeyed,  and  (iii)  the  form  of  the 
equations  remains  the  same  as  the  conventional  transfer  equations  so  that 
numerical  solutions  are  calculated  in  the  same  manner.  In  addition,  the 
dense  media  transfer  equation  takes  into  account  (i)  scattering  by 
correlated  scatterers,  (ii)  pair-distribution  function  of  scatterer 
positions,  and  (iii)  the  effective  propagation  constant  of  a  dense  medium. 

Numerical  solutions  of  the  dense-medi3  radiative  wave  equations  are 
illustrated  as  a  function  of  incident  angles,  scattered  angles,  and 
physical  parameters  of  the  medium.  We  have  derived  the  radiative  wave 
equations  for  the  scalar  wave  problem,  for  the  vector  electromagnetic 
propagation  of  the  active  remote  sensing  problem,  and  for  the  passive 
remote  sensing  problem.  The  results  of  passive  remote  sensing  have  found 
good  agreement  with  measurements  of  snow. 

Polarimetric  signatures  are  also  calculated.  Presently,  we  are 
extending  the  model  to  study  high  frequency  effects  as  well  as  media  with 
multi-species  of  particles. 
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B.  Backscattering  Enhancement 

Backacattering  enhancement  has  been  observed  in  laboratory 
controlled  experiments  of  scattering  from  discrete  particles.  The  effect 
is  not  contained  in  the  ladder  approximation  nor  in  the  radiative  transfer 
theory  because  the  cyclical  diagrams  which  contribute  significantly  in  the 
backscattering  direction  have  been  ignored.  We  have  first  examined 
backscattering  enhancement  by  a  second- order  cyclical  theory  and  have 
completed  a  study  based  on  the  resummation  of  all  the  cyclical  diagrams. 
For  the  case  of  small  scatterers  or  Rayleigh  scatterers,  the  theory 
predicts  a  sharp  peak  in  the  backscattering  direction  of  angular  width  of 
the  order  2K"/K’  where  K*  and  K"  are  the  real  and  imaginary  parts, 
respectively,  of  the  effective  propagation  constant.  A  modified  cyclical 
resummation,  which  includes  the  correlation  of  particle  positions  for  dense 
media  with  small  particles,  has  also  been  performed. 

For  the  case  of  moderate-size  particles,  the  scattering  is 
non- isotropic.  We  studied  the  backscattering  enhancement  of  waves  by 

nonisotropic  scatterers.  Multiple  scattering  effects  are  included  by 
examining  the  summation  of  all  the  ladder  terms  and  all  the  cyclical  terms. 
If  the  observation  angle  is  in  the  neighborhood  of  the  backscattering 
direction,  it  is  shown  that  both  summations  can  be  related  to  the 
unidirectional  point  source  Green's  function  of  the  transport  equation. 
For  the  case  of  small  albedo  or  small  optical  thickness,  the  second-order 
theory  is  applied  to  calculate  the  Green's  function.  The  angular  width  of 
backscattering  enhancement  in  this  ca3e  is  of  the  order  of  the  coherent 
wave  attenuation  rate  divided  by  the  wavenumber.  For  the  case  of  large 
albedo  and  large  optical  thickness,  the  diffusion  approximation  is  used  to 
calculate  the  Green's  function.  For  this  case,  the  angular  width  is  of  the 
order  of  the  transport  rate  divided  by  the  wavenumber.  The  transport  rate 
is  equal  to  the  product  of  the  coherent  wave  attenuation  rate  and  one  minus 
the  mean  cosine  of  the  scattering  angle.  Hence,  the  angular  width  is 
substantially  smaller  for  particles  with  dominant  forward  scattering  and  is 
in  good  agreement  with  experimental  observations.  Presently,  we  are 
extending  the  diffusion  approximation  to  the  case  of  vector  cyclical 
equations  and  also  that  of  beam  propagation.  The  effects  of 
pair-distr ibution  functions  are  also  investigated. 


C.  Inverse  Problems  in  Random  Media 

The  inverse  source  problem  of  the  scalar  wave  equation  for  a 
monochromatic  source  is  generalized  to  the  case  of  an  inhomogeneous 
attenuative  medium.  The  attenuative  medium  can  be  a  lossless  random  medium 
or  a  lossy  deterministic  medium  in  which  the  coherent  field  attenuates. 
Two  generalized  holographic  imaging  equations  are  obtained  based  on  the 
usage  of  two  Green'3  function  G+*  and  G_.  The  kernel  of  the  first  equation 
with  G  *  is  Hermitian  and  depends  on  the  location  and  shape  of  the 
recording  surface  while  the  equation  G-  has  a  non-Hermitian  kernel  that  is 
independent  of  the  recording  surface.  The  solutions  of  the  integral 
equations  are  investigated.  The  nonuniqueness  of  the  solutions  are  also 
related  to  the  nonradiating  sources  and  to  the  minimum  energy  solution. 
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D.  Conventional  Transport  Theory 


Me  have  investigated  the  small-angle  approximation  to  discrete 
random  media  for  the  transport  theory.  Studies  were  performed  for  both 
plane  wave  and  beam  wave  solutions.  It  is  found  that  for  particles  with  an 
index  of  refraction  close  to  1  and  a  size  parameter  larger  than  10,  the 
agreement  is  excellent.  However,  for  optical  propagation  through  fog,  the 
difference  between  the  radiative  transfer  and  small-angle  approximation  can 
be  more  than  30%. 
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Radiative  Wave  Equations  for  Vector  Electromagnetic  Propagation 
in  Dense  Nonienuous  Media 

Leung  Tseng  and  Akira  Ishimam 

Department  of  Electrical  Engineering 
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Abstract-  A  set  of  radiative  wave  equations  including  all  four  Stokes  parameters  is  derived 
for  vector  electromagnetic  wave  propagation  in  dense  nontenuous  media.  The  derivation  is 
based  on  the  quasi-crystalline  approximation  with  coherent  potential  on  the  first  moment 
of  the  field,  and  the  modified  ladder  approximation  on  the  second  moment  of  the  held. 
These  two  approximations  are  shown  to  be  energetically  consistent  for  dense  nontenuous 
media.  To  simplify  the  derivation  of  the  radiative  wave  equations,  the  model  of  small 
spherical  scatterers  is  used.  The  derived  radiative  wave  equations  assume  the  same  form 
as  the  classical  radiative  transfer  equations.  However,  the  relations  of  the  extinction  rate, 
the  albedo  and  the  phase  matrix  to  the  physical  parameters  of  the  media  include  the 
effects  of  dense  media  and  can  be  different  from  the  classical  relations  of  independent 
scattering. 


INTRODUCTION 

In  a  dense  medium,  the  particles  occupy  an  appreciable  fractional  volume.  In 
a  nontenuous  medium,  the  dielectric  properties  of  the  particles  are  significantly 
different  from  that  of  the  background  medium.  In  recent  years,  studies  have 
been  made  in  the  propagation  and  scattering  of  waves  in  dense  nontenuous  media 
with  applications  to  geophysical  terrain  such  as  snow,  ice-covered  land,  soil  and 
vegetation  [1-9].  In  a  dense  nontenuous  medium,  the  assumption  of  independent 
scattering,  that  is  often  used  in  radiative  transfer  theory  [10],  is  not  valid.  The 
model  of  a  continuous  random  medium,  that  has  been  applied  to  turbulent  media 
with  a  small  refractive  index  fluctuation  [111  is  not  appropriate  and  the  discrete 
scatterer  model  is  to  be  preferred.  Differences  between  continuous  medium  theory 
and  discrete  scatterer  theory  have  also  been  discussed  [12] . 

To  reconcile  the  radiative  transfer  theory  and  the  analytic  wave  theory,  it  is 
possible  to  derive  radiative  transfer  type  equations  from  the  second  moment  equa¬ 
tion  of  the  analytic  wave  theory  [8,12].  These  types  of  equations  are  similar  in 
form  to  the  classical  transfer  equation.  We  call  them  radiative  wave  equations  to 
distinguish  them  from  the  classical  radiative  transfer  equations  because  (1)  they 
are  derived  from  wave  theory  and  not  heuristic  and  (2)  they  often  correct  the 
deficiencies  of  the  classical  transfer  equation  when  the  latter  is  not  valid.  In  such 
cases,  major  differences  exist  between  the  radiative  wave  equations  and  the  das- 
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Passive  Remote  Sensing  of  Dense  Nontennous  Media 

L.  Tsang 
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Abstract'  The  thermal  microwave  —"«“*<»"  of  madia  with  dense  distribution  of  particles 
is  studied.  A  general  relation,  that  holds  for  media  with  nonuniform  at  well  as  uniform 
temperature  distribution,  is  first  derived  between  active  and  passive  remote  sensing.  The 
brightness  temperature  of  passive  sensing  is  proportional  to  the  integration  of  the  product 
of  the  temperature  distribution  and  the  divergence  of  the  Pointing's  vector  of  active  sens¬ 
ing.  The  governing  radiative  transfer  equations  for  passive  remote  tensing  of  dense  media 
are  then  derived  by  reference  to  the  corresponding  radiative  wave  equations  for  active 
sensing.  The  radiative  wave  equations  are  based  on  quasicrystalline  approximation  with 
coherent  potential  for  the  first  moment  of  the  field  and  a  modified  ladder  approximation 
for  the  second  moment  of  the  field.  The  numerical  results  for  classical  radiative  transfer 
theory  and  the  dense  media  transfer  theory  are  compared.  The  theory  is  also  used  to 
compare  with  brightness  temperature  measurements  over  a  snow  field. 

I.  INTRODUCTION 

In  recent  years,  studies  have  been  made  in  the  propagation  and  scattering  of 
waves  in  dense  nontenuous  media  with  applications  to  geophysical  terrain  such  as 
snow,  ice-covered  land,  soil  and  vegetation  [1-10].  In  a  dense  nontenuous  medium, 
the  assumption  of  independent  scattering,  that  is  often  used  in  classical  radiative 
transfer  theory,  [11)  is  not  valid.  To  reconcile  the  radiative  transfer  theory  and  the 
analytic  wave  theory,  it  is  possible  to  derive  radiative  transfer  type  equations  from 
the  second  moment  equation  of  the  analytic  wave  theory  [9,  12].  We  call  them 
radiative  wave  equations  to  distinguish  them  from  the  classical  radiative  transfer 
equations  because  (1)  they  are  derived  from  wave  theory  and  not  heuristically, 
and  (2)  they  often  correct  the  deficiencies  of  the  classical  transfer  equation  when 
the  latter  is  not  valid. 

In  this  paper,  a  set  of  radiative  wave  equations  is  derived  for  passive  remote 
sensing  of  media  with  a  dense  distribution  of  particles.  A  general  relation  is 
first  derived  between  active  and  passive  remote  sensing,  that  holds  for  media  with 
nonuniform  as  well  as  uniform  temperature  distribution.  By  using  the  fluctuation- 
dissipation  theorem  [9]  and  the  symmetry  relation  of  Green’s  functions,  it  is  shown 
that  the  brightness  temperature  of  passive  sensing  is  proportional  to  the  integra¬ 
tion  of  the  product  of  the  temperature  distribution  and  the  divergence  of  the 
Poynting’s  vector  of  active  sensing.  The  governing  radiative  wave  equations  for 
passive  sensing  are  next  derived  by  reference  to  the  corresponding  equations  that 
were  derived  for  active  sensing  in  a  previous  paper  [13].  The  radiative  wave  equa- 
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The  inverse  source  problem  of  the  scalar  wave  equation  for  a  monochromatic  source  is  generalized  to  the  case  of  an 
inhomogeneous  attenuative  medium.  The  attenuative  medium  can  be  a  loaay  deterministic  medium  or  a  lossless 
random  medium  in  which  the  coherent  field  attenuates.  Two  generalized  holographic  imaging  aquations  are 
obtained  that  are  based  on  the  use  of  two  Green's  functions,  G  ♦*  and  G-.  The  kernel  for  the  first  aquation  with  G+* 
is  Hermitian  and  depends  on  the  location  and  the  shape  of  the  recording  surface,  whereas  the  equation  with  G-  has  a 
non- Hermitian  kernel  that  is  independent  of  the  recording  surface.  The  solutions  of  the  integral  equations  are 
investigated.  The  nonuniqueneee  of  the  solutions  are  also  related  to  the  nonradiating  sources  and  to  the  minimum 
energy  solution. 


INTRODUCTION 


PROBLEM  STATEMENT 


A  problem  that  arises  in  optics  and  acoustics  is  the  inverse 
source  problem,1*4  the  solution  of  which  consists  of  deducing 
the  source  from  measurements  of  the  radiated  field  per¬ 
formed  outside  the  volume  containing  the  source.  The  in¬ 
verse  source  problem  is  usually  formulated  in  terms  of  inte¬ 
gral  equations  relating  the  source  to  some  quantity  calculat¬ 
ed  from  the  measured  fields.  Such  equations  are  sometimes 
called  generalized  holographic-imaging  equations.  Studies 
have  also  been  made  on  the  uniqueness  of  the  solution  of  the 
integral  equation,  which  has  been  related  to  the  possible 
existence  of  nonradiating  sources.3*6  Previous  treatment  of 
the  inverse  source  problem  has  assumed  the  medium  to  be 
.  imogeneous  and  lossless.  Recently,  an  extension  has  been 
made  to  the  case  in  which  the  medium  is  inhomogeneous  and 
nonabsorbing.' 

In  many  practical  situations,  the  background  medium  is 
lossy.  Absorption  effects  are  often  important  in  optics,  mi¬ 
crowaves.  and  acoustics.  For  the  case  of  wave  propagation 
in  random  media.39  the  coherent  wave  propagates  in  an 
effective  medium  and  attenuates  with  distance.  In  this  pa¬ 
per  we  generalize  the  results  of  the  inverse  source  problem  of 
Refs.  4  and  7  to  an  inhomogeneous  attenuative  medium. 
Two  integral  equations  are  obtained  that  are  baaed  on  the 
use  of  two  Green's  functions.  In  the  first  case,  the  backprop- 
agated  wave  attenuates  back  to  the  source  region.  The 
kernel  is  complex  Hermitian.  However,  it  is  dependent  on 
the  location  and  the  shape  of  the  recording  surface.  The 
mathematical  development  in  this  case  is  similar  to  that  of 
Ref.  7.  In  the  second  case  the  backpropagated  wave  grows 
back  to  the  source  region.  The  kernel  is  independent  of  the 
location  and  the  shape  of  the  recording  surface,  which  has  to 
lie  outside  the  source  region.  However,  the  kernel  is  com¬ 
plex  non-Hermitian.  The  solutions  of  the  two  integral 
equations  are  investigated.  The  nonuniqueness  of  the  solu¬ 
tions  is  also  related  to  the  nonradiating  sources  and  to  the 
minimum  energy  sources. 


We  consider  the  inverse  source  problem  for  a  localized 
monochromatic  source  p(r)  embedded  in  an  inhomogeneous 
attenuative  medium.  The  wave  equation  is 

[V2  +  *02n2(r)]iMr)  -  -/>(r),  (1) 

where  k0  is  the  free-space  wave  number.  The  index  of  re¬ 
fraction  is  a  complex  quantity  and  is  a  function  of  position: 

n(r)  -  n'(r)  +  in'(r),  (2) 

with  n'(r)  >  0  and  n”(r)  >  0.  Superscript  primes  and  dou¬ 
ble  primes  are  used  to  denote  the  real  and  imaginary  parts, 
respectively.  The  source  region  is  r,  and  the  recording  sur¬ 
face  Z  is  outside  the  source  region.  The  wave  function  and 
its  normal  derivative  are  known  on  the  recording  surface. 

There  are  two  Green's  functions  G*(r,  rD  and  G-(r,  r') 
associated  with  the  medium.  The  Green's  function  G+  de¬ 
notes  an  outgoing  attenuative  wave  from  a  point  source,  and 
the  Green's  function  G-  represents  an  incoming  growing 
wave: 

[7"J  +  fc,,;rt:ir>)G jr.  r'l  =  — 6i r  -  r'l.  131 

[V2  +  fe0-n-(r))G_(r.  r'l  =  -5<r  -  r'l.  (4) 

It  then  follows  that  the  wave  function  ^  is 

i£(r)  »  I  dr'G.lr,  r'lplrl.  (5) 


For  the  case  in  which  the  medium  is  homogeneous  with  n(r> 
■  n0,  the  Green’s  functions  Go*  and  Go-  are,  respectively. 

expiiklr  -  rl) 

4*lr  -  rl 


Go+(r.  r'l 


(61 


G0_(r,  r'l 


expi-ifrlr  -  rl) 
4wir  -  r'l 


(71 
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ABSTRACT  —  The  phase  pe mutation  technique  is  a  method  for  treating  wave  scattering  from 
rough  surfaces.  Winebrenner  and  Ishimaro  showed  that  the  phase  perturbation  expressions  for 
die  reflection  and  tackscattering  coefficients  reduce  to  those  of  classical  field  perturbation  theory 
for  surface  roughness.  Their  analysis  also  indicated  that  die  reflection  coefficient  reduces 
exactly  and  the  tackscattering  coefficient  approximately  to  those  of  the  Kirchhoff  approximation 
for  gently  undulating  surfaces.  We  numerically  verify  these  results  for  a  one-dimensional  surface 
with  Gaussian  surface  spectrum  satisfying  Dirichlet  boundary  conditions.  We  find  that  the  phase 
perturbation  backscattering  coefficient  reduces  numerically  to  that  of  Kirchhoff  over  a  large  range 
of  parameter  values  when  the  Kirchhoff  approximation  is  considered  to  be  valid.  In  addition,  the 
phase  perturbation  results  differ  from  die  field  perturbation  and  Kirchhoff  results  when  neither 
classical  method  is  expected  to  give  accurate  results.  Energy  conservation  is  also  considered. 

perturbation  methods 
Kirchhoff  approximation 
rough  surface  scattering 
wave  scattering 


random  Gaussian  surfaces 


BACKSCATTERING  ENHANCEMENT  OF  RANDOM  DISCRETE 
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Abstract: 

The  scattering  of  waves  by  particles  of  moderate  size  is 
nonisotropic.  In  this  paper,  the  backscattering  enhancement  of 
scattering  of  waves  by  nonisotropic  scatterers  is  studied.  Multiple 
scattering  effects  are  included  by  examining  the  summation  of  all  the 
ladder  terms  and  all  the  cyclical  terms.  If  the  observation  angle  is  in 
the  neighborhood  of  the  backscattering  direction,  it  is  shown  that  both 
summations  can  be  related  to  the  unidirectional  point  source  Green's 
function  of  the  transport  equation.  For  the  case  of  small  albedo  or 
I  small  optical  thickness,  the  second  order  theory  is  applied  to  calculate 

the  Green's  function.  The  angular  width  of  backscattering  enhancement  in 
this  case  is  of  the  order  of  the  coherent  wave  attenuation  rate  divided 
by  the  wavenumber.  For  the  case  of  large  albedo  and  large  optical 
thickness,  the  diffusion  approximation  is  used  to  calculate  the  Green's 
function.  For  this  case,  the  angular  width  is  of  the  order  of  the 
transport  rate  divided  by  the  wavenumber.  The  transport  rate  is  equal  to 
the  product  of  the  coherent  wave  attenuation  rate  and  one  minus  the  mean 
cosine  of  the  scattering  angle.  Hence,  the  angular  width  is  substan¬ 
tially  smaller  for  particles  with  dominant  forward  scattering  and  is  in 
good  agreement  with  experimental  observations. 
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BI1/P4-5  RADI ATI VI  TRANSFER  EQUATIONS  MODIFIED  FOR 
1020  APPLICATION  TO  DENSE  NONTENUOUS  MEDIA 

Leung  Tung  and  Akira  Ishlmaru 
Department  of  Elaetrleal  Englnaarlng 
Unlvaraity  of  Washington 
Saattla,  WA  98195 

Tha  classical  radlatlva  transfer  aquation  has  bean  used 
extensively  In  studying  volume  scattaring  in  raaota  sensing 
problaas.  Tha  advantages  of  using  the  radlatlva  transfer 
theory  are  that  (1)  multiple  scattering  of  the  Incoherent 
Intensities  are  Included  In  the  theory,  (11)  energy  conser¬ 
vation  Is  satisfied,  and  (111)  nuaericel  solutions  are 
tractable.  However,  the  transfer  equation,  in  Its  classical 
fora.  Is  not  applicable  for  dense  nontenuous  media. 

In  a  dense  medium,  the  particles  occupy  an  appreciable 
fractional  volume.  In  a  nontenuous  medium,  the  permittivity 
of  the  particles  is  significantly  different  from  that  of  the 
background  medium.  In  a  dense  nontenuous  medium,  the 
assumption  of  Independent  scattering  of  the  particles  no 
longer  holds  so  that  the  classical  transfer  equation  Is  not 
valid.  In  this  paper,  we  derive  a  set  of  "new"  transfer 
equations  from  analytic  wave  theory.  The  approximations 
made  on  the  analytic  wave  theory  are  (1)  quasl-crystalllne 
approximation  with  coherent  potential  on  the  first  moment  of 
the  field  and  (11)  a  modified  ladder  approximation  for  the 
second  moment  of  the  field.  These  two  approximations  are 
shown  to  satisfy  energy  conservation.  The  differences 
between  the  "new"  set  of  transfer  equations  and  the 
classical  ones  will  be  discussed.  Both  the  scalar  wave 
equations  and  the  vector  electromagnetic  propagation  will  be 
studied. 


PASSIVE  REMOTE  SENSING  OF  TERRAIN  WITH  A 
DENSE  DISTRIBUTION  OF  PARTICLES 

Leung  Tseng 

Department  of  Elaotrloal  Englnaarlng 
University  of  Washington 
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Passlva  remote  sensing  has  been  applied  to  terrain  media 
where  particles  are  densely  distributed.  The  radiative  transfer 
equation,  in  its  elassloal  form,  has  been  commonly  used  to  study 
the  volume  scattering  effects  of  these  problems.  However,  it  has 
been  demonstrated,  both  experimentally  and  theoretically,  that 
Independent  scattering  is  not  true  in  media  with  a  dense 
distribution  of  particles,  and  the  classical  radiative  transfer 
equation  is  not  valid. 

Recently,  a  set  of  "new"  transfer  equations  has  been  derived 
from  analytic  wave  theory  for  dense  media.  The  differences 
between  the  "new"  set  of  transfer  equations  and  the  classical 
ones  will  be  discussed.  The  new  equations  still  preserve  the 
advantages  that  (i)  multiple  scattering  of  the  Incoherent 
intensities  are  included,  (ii)  energy  conservation  is  satisfied, 
and  (lii)  numerical  solutions  are  tractable.  Numerical  results 
of  the  brightness  temperatures  based  on  the  theory  will  be 
illustrated  using  the  physical  parameters  of  the  terrain  media. 
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AimMpnmc  piniciM  Anver  w  graon  ctimtti 
by  absorbing  and  scattering  atactromagnaUc  radla- 
non.  isnowwogi  or  rot  optical  propvrovv  at  so 
lar  wavelength*  (0.3-4  Mm)  and  tarraatrlal  w ava- 
langBi*  (8-20  Mm)  la  naadad.  Photoacouatlc  mea¬ 
surement*  uaing  broadband  vWbla  light  have  baan 
pravWualy  mad*.  raaultlng  in  vahiaa  tor  lha  atmo¬ 
spheric  aaroaol  absorption  coefficient  of  1  X  10~’ 
m'V  Tha  photoacoustic  instrumant  has  baan 
modified  for  narrowband  maasuramants  in  lha  ul- 
traviolat  through  inlrarad  rag  Ions  Absorption  co- 
sfficiant  spectra  hava  baan  obtainad  for  atmo¬ 
spheric  particulates  collected  at  mountain  altitudes 
and  in  Tucson.  Arizona.  This  information  has 
been  used  for  identification  of  the  absorbing  mala¬ 
rial. 

Spectral  differences  are  found  for  samples  col¬ 
lected  in  urban  and  remote  locations.  Differences 
are  aiso.seen  m  spectra  for  small  particles  (<0.4 
Mm)  aiM  large  particles  (>0.4  pm).  ( 12  min) 

1  y  C.  Ramsey -Bell  and  Q.  Couture,  Appl.  Opt. 
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-/FL3  Comparisons  between  tha  small  angle  ap¬ 
proximation  and  tha  numerical  solution  for  radte- 
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sometmee  been  used  for  dMoiete  random  medto 
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Laser  beam  propagation  trou^i  fog  or  doud  Is 
a  problem  of  considerable  Interest  Aitiou^inu- 
marieel  sotuUone  of  tie  radtatve  tanalar  aquation 
in  tvaMabto  lor  Vo  com  of  pioo  wowo  Inoldwioo, 
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The  smell-angle  approxlmetlon  for  tie  radiative 
transfer  tieory  is  based  on  tha  assumption  that  tie 
scattered  wave  la  confined  within  a  smalt  angla  In 
tie  forward  dfreetton.  The  assumption  la  general- 
ly  valid  for  the  turbulence  case  In  which  tw  hrbu- 
lenoe  size  la  much  greeter  tiew  tie  wevetengti  and 
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The  fundamental  solution  la  a  result  of  savar« 
different  analytical  methods.  To  date,  only  a  com. 
parison  of  tie  edtWHetlon  index  and  tie  spactiaw 
of  intensity  fluctuations  (or  aqurvelantty  the  cone, 
latton  function)  has  baan  made,  wtti  numerical 
simulation  and  experimental  data-  We  extend  tha 
comparison  of  the  two  scale  results  to  the  Ml 
spatial  behavior  of  tha  fourth  moment.  Wacom- 
pared  our  results  wtti  nunaricel  sotuUone  obtained 
by  tie  fimtedHfarence  method  and  found  ttiat  tw 
typical  scales  of  the  fotrti  moment's  evolution  are 


Properties  such  as  symmsty/asymmeby  and 
dacorralaUon  are  also  clearly  dNpteyed.  (I2mla) 
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for  the  beam  wave  case  because  of  tha  formidabts 
amount  of  computation  Uma  required.  Tha  small 
angle  approxlmetlon  has  baan  shown  to  be  useful 
for  lightwave  propagation  In  babulam  motto  where 
the  variance  of  Index  of  refraction  fluctuations  la 
small.  However,  its  applicability  to  Wbid  meets 
with  discrete  scattarers  la  not  dear  because  of  tie 
much  larger  contrast  in  die  Max  of  refraction. 
Tha  small-angle  approximation  Is  uaad  to  calculate 
the  beam  wave  solution  throu0t  foggy  atmos¬ 
pheres  at  visible  wavelengths.  The  Hsnyey- 
Qreenstein  phase  function  Is  used  ®  describe  the 
scattering  characteristics  of  Use  panic  las.  The 
mutual  coherence  function  is  first  calculated  and 
than  the  Intensity  and  His  specific  Intensities  are 
calculated.  Numerical  results  are  illustrated  as  a 
function  of  optical  thickness  and  Hie  mean  cosine 
of  scattering  angles.  They  are  compared  with 
available  exact  numerical  solutions  of  tha  trans¬ 
port  equations  and  are  also  used  to  explain  Hie 
results  observed  In  laboratory  experimental  data. 

(12  min) 


FL5  Power  spectra  of  turbulence- Induced 
wavefront  aberrations  and  optical  path  difference 

R.  R.  BUTTS.  J.  F.  RIKER.  RSD  Associates.  105  E 
Vermijo.  Colorado  Springs.  CO  80903. 

Whan  a  point  source  radiates  wavefronts  which 
traverse  a  turbulent  atmosphere.  Hie  received 
wavefronts  will  contain  aberrations  which  vary  in 
time.  In  addition,  if  sevaral  point  sources  radiate 
in  an  array  conflgiratton,  Hie  resulting  wavefronts 
will  contain  optical  path  differences.  As  a  result. 
Hie  phase  profile*  incident  on  an  unobscured  col¬ 
lecting  apartura  will  b*  distorted  in  a  time-varying 
pattam. 

Ws  hava  expanded  Hi*  phase  profile  in  a  basis 
sat  of  Zamik*  polynomial*  end  computed  Hi*  pow¬ 
er  spectrum  appropriate  to  each  polynomial  coeffi¬ 
cient.  The  technique  used  incorporates  Parsa- 
val's  identity  to  allow  evaluation  of  certain  dinar 
products  In  Hi*  frequency  domain.  TW*  allow*  us 
to  simplify  Hi*  analysis  considerably  and  to  extend 
it  to  arbitrary  orders  of  basis  sat  polynomial*. 

W*  have  written  a  computer  program  to  quickly 
•viiuitw  9m  ipicti  for  rvtotivvty  towordif  ttoor- 
r Attorn  in  wtom  conditions,  snd  wo  prooont  somo 
of  thow  route.  (12  min) 


Tho  fUM  angular  dMbudont  for  mooring  of  • 
mm  mu  ■  ix  i  mil  pnnt  nw  oy  i  pons  coy  oo^ 
ducting  Aphvrv,  obtoinod  by  nrarted  mtimAbr 
of  tho  oxtet  poVil  wvw  ixpindon,  w  oomporod 
with  tho  prodtedons  of  oompltx  sngulor  momsv 
turn  (CAM)  theory-  Tho  olio  poromotor  x  ■  fca 
(whoro  k  -  wovonumbor.  «  •  tphoro  radtat) 
rongot  bom  largo  vatuoa  down  to  vahiaa  of  ordsr 

unnj,  w  oMnw  hm  mTwB  or  or 

CAM  theory.  A  complete  eat  of  scattering  data, 
I.*.,  Hi*  polarized  ditanaltlaa  and  Hi*  phaae  defer¬ 
ence  between  the  two  pofarizatlona.  la  plottad. 

In  the  panunbra  region,  beyond  Hie  forward 
diffraction  peak  (Airy  pattam)  and  up  to  scattering 
angla*  9  at  order  CAM  theory  leads  to  a  new 

and  improved  version  of  Fock's  theory  of  diffrac¬ 
tion.  Hi*  uniform  Fock  approximation.  Beyond 
Hill  region,  scattaring  la  dominated  by  Hie  WKB 
approximation,  but  surfaca-wav*  contributions 
also  play  an  important  role,  especially  In  the  tower 
range  of  value*  of  x  and  9. 

The  asymptotic  predictions  of  CAM  theory  ar*  In 
excellent  agreement  with  Hie  exact  angular  scat¬ 
tering  panems.  ( 12  mki) 

FL9  Lota  of  Image  resolution  In  haz*  and  duet 

P  L  WALKER.  U  S.  Naval  Weapons  Center.  China 
Lake.  CA  93555 

Scattering  ol  light  by  particulate  maner  in  the 
atmosphere  generates  a  glow  or  aureol  around 
pointlike  sources,  such  as  street  lights  seen  from  a 
distance.  A  scene  can  be  thought  ol  as  an  ensem¬ 
ble  of  pointlike  sources  or  pixels,  each  having  an 
aureol  associated  with  it.  This  phenomena 
causes  images  of  Hie  scene  fo  be  smeared:  hence, 
the  presence  of  Hie  scetterers  in  the  path  of  propa¬ 
gation  causes  loss  of  image  resolution  a  compu¬ 
tational  method  we*  previously  reported’  whereby 
a  Monte-Carlo  code  coupled  with  an  optica  design 
program  waa  used  to  generate  atreois  at  Hi*  kn¬ 
ag*  plan*  of  a  detector  This  procedure  now 
include*  a  Gaussian  MTF  for  Hi*  camera,  so  that 
the  optica  of  a  camera  operating  in  a  scattaring 
environment  can  b*  completely  simulated. 

Auraot*  generated  by  Hit*  method  ar*  compared 
with  tho**  by  an  analytical  method  propoaad  by 
Lutomlrskl.,  It  Is  determined  In  what  condkton*  of 
optical  depth  and  paakadnaaa  of  phase  function 
HiatwomaHiedsgfvaHiasamarssuNB.  ins  1991 
paper  Kopatra  at  at*  observed  enhanced  propaga¬ 
tion  wtti  increasing  asvstsngM  In  Hi*  near  to  mid 
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BIO-4  3ACKSCATTERING  ENHANCEMENT  FOR  RANDOM  DISCRETE 

0940  SCATTERERS  OF  MODERATE  SIZE 

Leung  Tsang  and  Akira  Ishimaru 
Department  of  Elactrical  Engineering 
University  of  Washington 
Seattle,  Washington  98195 

Backscattertng  enhancement  for  random  discrete 
scatterers  was  previously  studied  by  calculating  the 
contributions  of  cyclical  scattering  terns  in  multiple 
scattering  processes.  The  method  was  applied  to  the  case  of 
small  particles,  where  it  was  shown  that  the  cyclical 
scattering  processes  exhibit  a  peak  in  the  backscattering 
direction  with  angular  width  2K"/K’  with  K'  and  K"  being  the 
real  and  imaginary  parts  of  the  effective  propagation 
constant.  In  this  paper,  the  method  is  extended  to  the  case 
of  scatterers  comparable  to  or  larger  than  a  wavelength. 
The  problem  of  a  plane  wave  normally  Incident  on  a  half 
space  of  dielectric  spheres  is  considered.  The  problem  is 
of  particular  interest  because  besides  the  peak  due  to 
cyclical  scattering,  there  can  also  be  a  peak  due  to  single 
scattering  in  the  Mle  phase  function.  Numerical  results  are 
illustrated  as  a  function  of  particle  size  and 
concentration.  The  relative  importance  of  the  two  peaks 
will  be  discussed. 
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CONVENTIONAL,  AND  MEN  RADIATIVE  TRANSFER  EQUATIONS  FOR 
VOLUME  SCATTERING  IN  GCOFNISICAL  MEDIA 


laung  TMng  and  Bohang  Wan 
Dapartasnt  of  Elactrteal  Engineering 
University  of  Washington 
Seattle,  Washington  9*195 


Recently,  a  set  of  no*  radiative  tranafar  equations  has  been 
derived  fro*  wave  theory  for  electromagnetic  wave  propagation  In 
dansa  geophysical  aadla.  Tha  aotlvatton  la  that  conventional 
tranafar  aquatlona  assuaa  Lndapandant  >  oat  tar  me  wltlch  la  not  valid 
in  danaa  aadla.  Tha  dlffaraneas  Mtwaan  tha  eonvantlonal  and  naw 
tranafar  aquatlona  will  be  dlseusaad.  Tha  naw  aquatlona  still 
praaarva  tha  advantagaa  that  (1)  aultlpla  soattarlns  of  tha 
ineonarant  lntanaltlaa  a ra  includad,  (11)  energy  eonaarvatlon  and 
raeiproclty  ara  obayad,  and  (111)  tna  fora  of  tha  aquatlona  remains 
tha  saaa  as  tha  oonvantlonal  tranafar  aquatlona  so  tnat  numerical 
solutions  ara  calculatad  In  tha  saaa  aartnar.  In  addition,  tna 
danaa  aadla  tranafar.  aquation  takas  into  acoount  (l)  scattaring  by 
corralatad  scattarars,  (11)  pair-distribution  function  or  scattarar 
positions,  and  (111)  tha  affactlva  propagation  constant  of  a  dansa 
medium.  In  this  papar ,  wa  also  ir.ow  numerical  rasulta  baaad  on  tha 
naw  transport  aquatlona  using  -.na  physical  paraaatars  of  tna 
tarraln  aadla  and  coapara  than  »;-.n  that  of  tha  eonvantlonal  thaory 
and  axpanaantal  data  in  both  act.ve  and  paaaiva  raaota  sanslng. 
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Ahetruct 

We  iii  ii  eirigam  dw  polariatioa  progenies  o f  Km  scanned  field  according  to  a  discrete 
ifowr  modal.  Th»  modal  cm  apply  to  a  number  at  nuunl  madia:  in  particular.  we  use  it  to 
represent  flnt-  and  mold-year  sea  tea  at  L-  and  C-bend.  We  discuss  pondal  ways  to  infer 
geophysical  propernee  of  the  tcaaahag  medium  modal  using  compiec  information  about  the 
polarization  of  the  scattered  field  (e-g.,  position  of  polarization  nulls  and  degree  of  polarization), 
especially  in  caaae  wham  like-  and  croaa  polarized  bsdttcener  cross  sections  alone  do  not 
provide  much  geophysical  mfonnadon.  Although  dte  scattering  model  used  in  this  intial 
investigation  may  be  tooewhat  too  simple  to  predict  detailed,  quantitative  polarization 
propones  of  sea  ice.  it  does  provide  considerable  insight  and  guidance  in  choosing  poUnmetne 
quantities  to  examine  it  data  of  this  kind  become  available. 
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